Quantum pumping in closed systems, adiabatic transport, and the Kubo formula 
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Quantum pumping in closed systems is considered. We explain that the Kubo formula contains all 
the physically relevant ingredients for the calculation of the pumped charge (Q) within the frame- 
work of linear response theory. The relation to the common formulations of adiabatic transport 
and "geometric magnetism" is clarified. We distinguish between adiabatic and dissipative contri- 
butions to Q. On the one hand we observe that adiabatic pumping does not have to be quantized. 
On the other hand we define circumstances in which quantized adiabatic pumping holds as an ap- 
proximation. The deviation from exact quantization is related to the Thouless conductance. As an 
application we discuss the following examples: classical dissipative pumping by conductance control, 
classical adiabatic (non dissipative) pumping by translation, and quantum pumping in the double 
barrier model. In the latter context we analyze a 3 site lattice Hamiltonian, which represents the 
simplest pumping device. We remark on the connection with the popular S matrix formalism which 
has been used to calculate pumping in open systems. 



I. INTRODUCTION 

Linear response theory (LRT) [1-3] is the leading for- 
malism to deal with driven systems. Such systems are 
described by a Hamiltonian Ti{x) where x{t) is a set 
of time dependent classical parameters ("fields"). The 
Kubo formula is the corner stone of LRT. It allows the 
calculation of the response coefficients, and in particular 
the conductance matrix (G) of the system. If we know 
G, we can calculate the charge (Q) which is transported 
through the system during one cycle of a periodic driving. 
This is called pumping. 

Pumping of charge in mesoscopic [4] and molecular size 
devices is regarded as a major issue in the realization of 
future quantum circuits or quantum gates, possibly for 
the purpose of quantum computing. 



A. Model system 

In order to explain the motivation for the present work, 
and its relation to the published literature, we have to 
give a better definition of the problem. For presentation 
purpose we focus on a model system with a ring geome- 
try (Fig.l). The shape of the ring is controlled by some 
parameters xi and X2- These parameters can be gate 
voltages that determine the location of some boundaries, 
or the height of some barriers. The third parameter is 
the flux through the ring: 
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(1) 



We shall use units such that the elementary charge is 
e — 1. Note that the Hamiltonian Ti{xi{t) , X2{t) , X3{t)) 
has gauge invariance for cj) i—f (f> + 2-k. Another system 
with a ring topology is presents in Fig. lb, and its ab- 
straction is represented in Fig. 2c. The "dot" can be 
represented by an S matrix that depends on xi and X2- 
In Fig. Id also the flux x^ is regarded as a parameter of 
the dot. If we cut the wire in Fig. Id we get the open 
two lead geometry of Fig.le. Finally we can put many 



such units in series (no flux), hence getting the periodic 
system of Fig. If. In the latter case the Hamiltonian is 
invariant for unit translations, and therefore the quasi 
momentum ^ is a constant of motion. It follows that the 
mathematical treatment of a driven periodic structure 
reduces to an analysis of a driven ring system with flux. 



B. Classification of pumps 

"Pumping" means that net charge (or maybe better to 
say "net integrated probability current") is transported 
through the ring per cycle of a periodic driving. Using 
the common jargon of electrical engineering this can be 
described as AC-DC conversion. We distinguish between 

• pumping in open systems (such as in Fig. le). 

• pumping in closed systems (such as in Fig. Id) 

• pumping in periodic systems (such as in Fig. If) 

For a reason that was explained at the end of the previous 
subsection we regard the last category [5] as mathemat- 
ically equivalent to the second category. We also regard 
the first category [6-9] as a special (subtle) limit of the 
second category: in a follow up paper [10] we demon- 
strate that in the limit of open geometry the Kubo for- 
mula reduces to the 5-niatrix formula of Biittiker Pretre 
and Thomas [6]. 

There are works in the literature regarding "rectifi- 
cation" and "ratchets" [11]. These can be regarded as 
studies of pumping in periodic systems with the con- 
notation of having damped non-Hamiltonian dynamics. 
These type of systems are beyond the scope of the present 
Paper. There is also a recent interest in Hamiltonian 
Ratchets [12], which is again a synonym for pumping in 
periodic systems, but with the connotation of having a 
non-linear pumping mechanism. We are going to clar- 
ify what are the conditions for having a linear pumping 
mechanism. Only in case of linear pumping mechanism 
the Kubo formula can be used, which should be distin- 
guished from the non- linear mechanism of Ref.[12]. 



2 



C. Objectives 

The purpose of this Paper is to explain and demon- 
strate that the Kubo formula contains all the physically 
relevant ingredients for the calculation of the charge {Q) 
which is pumped during one cycle of a periodic driving. 
In the limit of a very slow time variation (small i), the 
emerging picture coincides with the adiabatic picture of 
Refs.[5, 13-15]. In this hmit the response of the system 
is commonly described as a non-dissipative "geometric 
magnetism" [15] effect, or as adiabatic transport. 

A major objective of this Paper is to bridge between 
the adiabatic picture and the more general LRT / Kubo 
picture, and to explain how dissipation emerges in the 
quantum mechanical treatment. For one-parameter driv- 
ing a unifying picture that bridges between the quan- 
tum mechanical adiabatic picture and LRT has been pre- 
sented in [16-19]. A previous attempt [20] had ended in 
some confusion regarding the identification of the linear 
response regime, while [15] had avoided the analysis of 
the mechanism that leads to dissipation in the quantum 
mechanical case. 

The presented (Kubo based) formulation of the pump- 
ing problem has few advantages: It is not restricted to 
the adiabatic regime; It allows a clear distinction between 
dissipative and adiabatic contributions to the pumping; 
The classical limit is manifest in the formulation; It gives 
a level by level understanding of the pumping process; It 
allows the consideration of any type of occupation (not 
necessarily Fermi occupation); It allows future incorpo- 
ration of external environmental influences such as that 
of noise; It regards the voltage over the pump as elec- 
tro motive force, rather than adopting the conceptually 
more complicated view [9] of having a chemical potential 
difference. 

Of particular interest is the possibility to realize a 
pumping cycle that transfers exactly one unit of charge 
per cycle. In open systems [7, 8] this "quantization" holds 
only approximately, and it has been argued [7] that the 
deviation from exact quantization is due to the dissipa- 
tive effect. Furthermore it has been claimed [7] that exact 
quantization would hold in the strict adiabatic limit, if 
the system were closed. In this Paper we would like to 
show that the correct picture is quite different. We shall 
demonstrate that the deviation from exact quantization 
is in fact of adiabatic nature. This deviation is related 
to the so-called "Thouless conductance" of the device. 



D. Examples 

We give several examples for the application of the 
Kubo formula to the calculation of the pumped charge Q: 

• classical dissipative pumping 

• classical adiabatic pumping (by translation) 

• quantum pumping in the double barrier model 



The last example is the main one. In the context of open 
geometry it is known as "pumping around a resonance" 
[8]. We explain that this is in fact an diabatic transfer 
scheme, and we analyze a particular version of this model 
which is represented by a 3 site lattice Hamiltonian. This 
is definitely the simplest pump circuit possible, and we 
believe that it can be realized as a molecular size device. 
It also can be regarded as an approximation for the closed 
geometry version of the two delta potential pump [8] . 



E. Outline 

In Section 2 we define the main object of the study, 
which is the conductance matrix G of Eq.(5). The con- 
ductance matrix can be written as the sum of a sym- 
metric (rj) and an anti-symmetric (B) matrices, which 
are later identified as the dissipative and the adiabatic 
contributions respectively. 

In the first part of the paper (Sections 2-8) we analyze 
the adiabatic equation (Section 3), and illuminate the 
distinction between its zero order solution (Section 4), 
its stationary first order solution (Section 5), and its non- 
stationary solution (Section 6). The outcome of the anal- 
ysis in Section 5 is Eq.(26) for the conductance matrix 
G. This expression is purely adiabatic, and does not give 
any dissipation. In order to get dissipation we have to 
look for a non-stationary solution. 

The standard textbook derivation of the Kubo formula 
(Eq.(29)) for the conductance matrix G implicitly as- 
sumes a non-stationary solution. We show how to get 
from it Eq.(30) for rj and Eq.(31) for B. The latter is 
shown to be identical with the adiabatic result (Eq.(26)). 
In Section 7 we further simplify the expression for rj lead- 
ing to the fluctuation-dissipation relation (Eq.(33)). 

The disadvantages of the standard textbook derivation 
of Kubo formula make it is essential to introduce a dif- 
ferent route toward Eq.(33) for rj. This route, which is 
discusssed in Section 8, explicitly distinguishes the dis- 
sipative cflcct from the adiabatic effcict. and allows to 
determine the conditions for the validity of either the adi- 
abatic picture or LRT. In particular it is explained that 
LRT is based, as strange as it sounds, on perturbation 
theory to inflnite order. 

In Section 9 we clarify the general scheme of the pump- 
ing calciflation (Eq.(39)). Section 10 and Section 11 give 
two simple classical examples. In Section 12 we turn 
to discuss quantum pumping, where the cycle is around 
a chain of degeneracies. The general discussion is fol- 
lowed by presentation of the double barrier model (Sec- 
tion 13). In order to get a quantitative estimate for the 
pumped charge we consider a 3 site lattice Hamiltonian 
(Section 14). 

The summary (Section 15) gives some larger perspec- 
tive on the subject, pointing out the relation to the S- 
matrix formalism, and to the Born-Oppenheimer picture. 
In the appendices we give some more details regarding 
the derivations, so as to have a self-contained presenta- 
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tion. 



II. THE CONDUCTANCE MATRIX 

Consider the Hamiltonian 'H{x{t)), where x{t) is 
a set of time dependent parameters ("fields"). For 
presentation, as weU as for practical reasons, we as- 
sume later a set of three time dependent parameters 
x{t) = {xi{t),X2{t),X3{t)). We define generalized forces 
in the conventional way as 



pK _ _ 



m 

dxk 



(2) 



Note that if xi is the location of a wall element, then is 
the force in the Newtonian sense. If X2 is an electric field, 
then is the polarization. If X3 is the magnetic field or 
the flux through a ring, then F^ is the magnetization or 
the current through the ring. 

In linear response theory (LRT) the response of the 
system is described by a causal response kernel, namely 



-< 



{t-t') Xj{t')dt' 



(3) 



where a''^ (t) = for t < 0. The Fourier transform 
of a'^^ (r) is the generalized susceptibility x'^H'^)- The 
conductance matrix is defined as: 



u^O U) Jo 



(4) 



Consequently, as explained further in Appendix B, the 
response in the "DC limit" {uj 0) can be written as 



(5) 



As an example for the applicability of this formula note 
the following standard examples for one parameter driv- 
ing: Let X = wall or piston displacement, then x = wall 
or piston velocity, G = friction coefficient, and F = —Gx 
is the friction force. Another standard example is x = 
magnetic flux, —x= electro motive force, G = electrical 
conductance, and hence F = —Gx is Ohm law. 

It is convenient to write the conductance matrix as 
gfcj = ry'^J + B^^ , where r}^^ = r]^^ is the symmetric part 
of the conductance matrix, while B'^^ = —B^'^ is the 
antisymmetric part. In case of having three parameters 
we can arrange the elements of the antisymmetric part as 
a vector B = (5^3^ ^3i^ ^12^ Consequently Eq.(5) can 
be written in abstract notation as 



{F) 



-rj ■ X — B Ads 



(6) 



where the dot product should be interpreted as matrix- 
vector multiplication, which involves summation over the 
index j. The wedge-product also can be regarded as a 



matrix-vector multiplication. It reduces to the more fa- 
miliar cross-product in case that we consider 3 param- 
eters. The dissipation, which is defined as the rate in 
which energy is absorbed into the system, is given by 

kj 

Only the symmetric part contributes to the the dissi- 
pation. The contribution of the antisymmetric part is 
identically zero. 



III. THE ADIABATIC EQUATION 

The adiabatic equation is conventionally obtained from 
the Schrodinger equation by expanding the wavefunction 
in the x-dependent adiabatic basis: 



= 5^a„(i) \n{x{t))) 



m j 

where following [13] we define 

/ d 



(8) 
(9) 

(10) 



(11) 



Differentiation by parts of dj {n{x)\in{x)) = leads to 
the conclusion that is a hermitian matrix. Note 

that the effect of gauge transformation is 



\n{x)) 



\n{x)) 

Aim + {djAn)6„ 



(12) 
(13) 



Note that the diagonal elements = A^^ are real, and 
transform as A^^ Ai + 9jA„. 

Associated with An{x) is the gauge invariant 2-form, 
which is defined as: 



= diAi - d.Al, 
= —2Jilm{din\djn) 



2 • 



(14) 
(15) 

(16) 



This can be written in an abstract notation as B = VAA. 
Using standard manipulations, namely via differentiation 
by parts of dj{n{x)\'H\m{x)) = 0, we get for n ^ m the 
expressions: 



ih 



R —R 



dU 



dxi 



R —R 



and hence 



B^- 



(17) 



(18) 
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IV. THE STRICTLY ADIABATIC SOLUTION, 
AND THE BERRY PHASE 

We define the perturbation matrix as 

Wnm = - ^ XjAf^^ for n ^ m (19) 



and = for n = rn. Then the adiabatic equation 

can be re- written as follows: 



dan 



UUy„ ' / T^ A \ ' TIT- 

— = --{En-xAn)an - J } ^ Wnma 



(20) 



If we neglect the perturbation W, then we get the strict 
adiabatic solution: 



-(/;E„(.(t'))dt'-/;;;'A„(.).d.) 



\n{xm (21) 



Due to An{x), we have the so called geometric phase. 
This can he gaiigc^d away unless we consider a closed 
cycle. For a closed cycle, the gauge invariant phase 
(l/h) ^ A ■ dx is called Berry phase. 

With the above zero-order solution we can obtain the 
following result: 



n{x) 

d 
dxk 



{n{x)\ n{x) \n{x)) (23) 



In case of the standard examples that were mentioned 
previously this corresponds to conservative force or to 
persistent current. From now on we ignore this trivial 
contribution to (F*) , and look for the a first order con- 
tribution. 



V. THE STATIONARY ADIABATIC 
SOLUTION: ADIABATIC TRANSPORT OR 
"GEOMETRIC MAGNETISM" 

For linear driving (unlike the case of a cycle) the An {x) 
field can be gauged away. Assuming further that the adi- 
abatic equation can be treated as parameter independent 
(that means disregarding the parametric dependence of 
En and W on x) one realizes that Eq.(20) possesses sta- 
tionary solutions. To first order these are: 



= \n)+ E 



En — E„ 



|m) 



(24) 



Note that in a fixed-basis representation the above sta- 
tionary solution is in fact time-dependent. Hence the no- 
tations \n{x{t))), \m{x{t))) and \tp{t)) are possibly more 
appropriate. 

With the above solution we can write {F'^) as a sum 
of zero order and first order contributions. From now on 



we ignore the zero order contribution, and go on with the 
first order contribution: 



Wrr 



En — E„ 



dH 



dxk 



m)+CC 



j \ m / 



(25) 



For a general stationary preparation, either pure or 
mixed, one obtains Eq.(5) with 



(26) 



where f{En) are weighting factors, with the normal- 
ization ^nf{Fn) = 1- For a pure state preparation 
/(£"„) distinguishes only one state n, while for canoni- 
cal preparation f{En) oc exp{—En/T), where T is the 
temperature. For a many-body system of non-interacting 
Fermions /(£"«) can be re-interpreted as the Fermi occu- 
pation function, so that .f{En) is the total number of 
particles. 

Thus we see that the assumption of a stationary first- 
order solution leads to a non-dissipative (antisymmet- 
ric) conductance matrix This is know as either Adiabatic 
Transport [5, 14] or "Geometric Magnetism" [13]. In the 
later sections we shall discuss the limitations of the above 
result. 



VI. THE NON-STATIONARY SOLUTION: THE 
KUBO FORMULA 

The Kubo formula is an expression for the linear re- 
sponse of a driven system that goes beyond the stationary 
adiabatic solution of the previous section. The Kubo for- 
mula has many type of derivations. One possibility is to 
use the same procedure as in Section 5 starting with 



n) 



E 



-iW„^ 



f 

Jo 



-iEmt 



m) 



For completeness we also give in Appendix A a simple 
version of the standard derivation, which is based on a 
conventional fixed-basis first-order-treatment of the per- 
turbation. The disadvantages are: 

• The standard derivation does not illuminate the un- 
derlaying physical mechanisms of the response. 

• The stationary adiabatic limit is not manifest. 

• The fluctuation-dissipation relation is vague. 

• The validity conditions of the derivation are not 
clear: no identification of the regimes. 
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For now wc go on with the conventional approach, but 
in a later section we refer to the more illuminating ap- 
proach of [20] and [16-19]. Then we clarify what is the 
regime (range of x) were the Kubo formula can be trusted 
[18, 19], and what is the sub-regime where the response 
can be described as non-dissipativc adiabatic transport. 

In order to express the Kubo formula one introduces 
the following definition: 



K'^{r) = -{[F\r),F\0)]) 



(27) 



We use the common interaction picture notation 
F'=(r) = e^wtffee-*"*, where H = n{x) with x = const. 
The expectation value assumes that the system is pre- 
pared in a stationary state (see previous section). It is 
also implicitly assumed that the result is not sensitive 
to the exact value of x. Note that K^^{t) has a well de- 
fined classical limit. Its Fourier transform will be denoted 

The expectation value {F^) is related to the driving 
x{t) by the causal response kernel a^^ [t — t'). The Kubo 
expression for this response kernel, as derived in Ap- 
pendix A, is 



a^^T) = e(r) K'^{t) 



(28) 



where the step function 6(t) cares for the upper cutoff of 
the integration in Eq. (3) . The Fourier transform of a^^ (r) 
is the generalized susceptibility x*-' (a;). The conductance 
matrix is defined as: 



poo 

= / K'\T)TdT (29) 
Jo 



This can be split into symmetric and anti-symmetric 
components (see derivation in Appendix C) as follows: 



^ij ^ i(G'J-G^') 



1 ^.^ lm[^ 

2 co-*0 U! 



■f 



2^ tj2 



(31) 



The antisymmetric part is identified [15] (Appendix C) 
as corresponding to the stationary solution Eq.(26) of the 
adiabatic equation. 



VII. THE EMERGENCE OF DISSIPATIVE 

RESPONSE, AND THE 
FLUCTUATION-DISSIPATION RELATION 

The Kubo formula for the symmetric part of the con- 
ductance matrix [t]^^) can be further simplified. If we 
take Eq.(30) literally, then 77*^ = due to the simple fact 
that we have finite spacing between energy levels (see 
[21] for a statistical point of view). But if we assume 
that the energy levels have some finite width F, then the 
smoothed version of K'^^{uj) should be considered. In 
common textbooks the introduction of F is "justified" by 



assuming some weak coupling to an environment, or by 
taking the limit of infinite volume. But we arc dealing 
with a strictly isolated finite system, and therefore the 
meaning of F requires serious consideration. We post- 
pone the discussion of this issue to the next section. 

If the smoothed version of K^^ (uj) should be used in 
Eq.(30), then it is possible to obtain 7]^^ from power 
spectrum C'*-' (w) of the fluctuations. This is called the 
Fluctuation-Dissipation relation. The spectral function 
C^^{u!) is defined as the Fourier transform of the sym- 
metrized correlation function 

C'^{t) = {1{F'{t)F'{0) + F'{0)F'{t))) (32) 

We use again the interaction picture, as in the definition 
of K'^^{t). Also this function has a well defined classical 
limit. 

There are several versions for the Fluctuation- 
Dissipation relation. The microcanonical version [15] has 
been derived using classical considerations, leading to 



1 1 



2 g{E) dE 



(33) 



In Appendix C wc introduce its quantum mechanical 
derivation. The subscript emphasizes that we assume a 
microcanonical state with energy E, and g{E) is the den- 
sity of states. The traditional version of the Fluctuation- 
Dissipation relation assumes a canonical state. It can be 
obtained by canonical averaging over the microcanonical 
version leading to: 



1 

2f 



0) 



(34) 



VIII. 



THE VALIDITY OF LINEAR RESPONSE 
THEORY AND BEYOND 



The standard derivation of the dissipative part of the 
Kubo formula, leading to the Fluctuation-Dissipation re- 
lation Eq.(33), is not very illuminating physically. More 
troubling is the realization that one cannot tell from the 

standard derivation what are the conditions for its va- 
lidity. An alternate, physically appealing derivation [16- 
19], is based on the observation that energy absorption 
is related to having diffusion in energy space [20]. In 
Appendix E wc outline the main ingredients of this ap- 
proach. 

It should be clear that the diffusion picture of Ap- 
pendix E holds only in case of chaotic systems. If this 
diffusion picture does not hold, then also the Kubo for- 
mula for T]^^ does not hold! Driven one-dimensional sys- 
tems are the obvious example for the failure of linear re- 
sponse theory (LRT). As in the case of the kicked rotator 
(standard map) [22] there is a complicated route to chaos 
and stochasticity: By increasing the driving amplitude 
the phase space structure is changed. If the amplitude 
is smaller than a threshold value, then the diffusion is 
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blocked by Kolmogorov-Arnold-Moser curves, and con- 
sequently there is not dissipation. Therefore the Kubo 
formula is not applicable in such cases. 

The following discussion of dissipative response as- 
sumes that we deal with a quantized chaotic system. We 
would like to discuss the reason and the consequences of 
having an energy scale T. In the standard derivation the 
assumption of having level broadening as if comes out of 
the blue. As we already noted it is customary in text- 
books to argue that either a continuum limit, or some 
small coupling to an environment is essential in order to 
provide T. But this is of course just a way to avoid con- 
frontation with the physical problem of having a driven 
isolated finite mesoscopic system. In fact the energy scale 
r is related to the rate (i) of the driving: 

r = i^^\x\y\A (35) 

where for simplicity of presentation we assume one pa- 
rameter driving. We use A to denote the mean level 
spacing, and a is the root mean square value of the ma- 
trix element Fnm between neighboring levels. In order 
to derive the above expression for F we have used the 
result of [19] (Sec. 17) for the "core width" at the break- 
time t = tprt of perturbation theory. The purpose of the 
present section is to give an optional " pedestrian deriva- 
tion" for r, and to discuss the physical consequences. 

Looking at the first order solution Eq.(24) of Section 5 
one realizes that it makes sense provided IWmnl A. 
This leads to the adiabaticity condition 

\i\ « ^ (36) 

If this condition is not satisfied one should go beyond 
first order perturbation theory (FOPT), in a sense to be 
explained below. Note that this adiabaticity condition 
can be written as F <C A. 

The adiabaticity condition Eq.(36) can be explained 
in a more illuminating way as follows: Let us assume 
that we prepare the system at time t = at the level 
|n). Using time dependent FOPT we find out that a 
stationary-like solution is reached after the Heisenbcrg 
time ts = 2TT?i/A. This is of course a valid description 
provided we do not have by then a breakdown of FOPT. 
The condition for this is easily found to be xt^ <C 5xc, 
where 6xc = A/a. This leads again to the adiabaticity 
condition Eq.(36). 

Another assumption in the derivation of Section 5 was 
that we can ignore the parametric dependence of £"„ and 
W on X. The adiabaticity condition t^ ^ Sxc/x mani- 
festly justify such an assumption: We should think of ta 
as the transient time for getting a stationary-like state, 
and we should regard 6xc as the parametric correlation 
scale. 

As strange as it sounds, in order to have dissipation, 

it is essential to have a breakdown of FOPT. In the lan- 
guage of perturbation theory this implies a required sum- 



mation of diagrams to infinite order, leading to an effec- 
tive broadening of the energy levels. By iterating FOPT, 
neglecting interference terms, we get a Markovian ap- 
proximation for the energy spreading process. This leads 
to the diffusion equation of Appendix E. This diffusion 
can be regarded as arising from Fermi-golden-rule tran- 
sitions between energy levels. A simple ad-hoc way to 
determine the energy level broadening is to introduce F 
as a lower cutoff in the energy distribution which is im- 
plied by Eq.(24): 

This constitutes a generalization of the well known proce- 
dure used by Wigner in order to obtain the local density 
of states [16, 17]. However, in the present context we 
do not get a Lorentzian. The width parameter F is de- 
termined self consistently from normalization, leading to 
Eq.(35) [disregarding numerical prefactor]. 

We can summarize the above reasoning by saying that 
there is a perturbative regime that includes an adiabatic 
(FOPT) sub-regime. Outside of the adiabatic sub-regime 
we need all orders of perturbation theory leading to 
Fermi-golden-rulc transitions, diffusion in energy space, 
and hence dissipation. Thus the dissipative part of Kubo 
formula emerges only in the regime F > A, which is just 
the opposite of the adiabaticity condition. The next ob- 
vious step is to determine the boundary of the perturba- 
tive regime. Following [18, 19] we argue that the rcriuircd 
condition is F <C A^. The bandwidth A^ oc ?i is defined 
as the energy width \En — Em\ were the matrix elements 
Fnm are not vanishingly small. If the condition F <C A;, is 
violated we find ourselves in the non-perturbative regime 
where the Kubo formula cannot be trusted [18, 19]. 

We still have to illuminate why we can get in the per- 
turbative regime a dissipative linear response in spite of 
the breakdown of FOPT. The reason is having a separa- 
tion of scales (A <C F <C A;,). The non-perturbative 
mixing on the small energy scale F does not affect 
the rate of first-order transitions between distant lev- 
els (F < \En - Em\ < Ab). Therefore Fermi golden rule 
picture applies to the description of the coarse grained 
energy spreading, and we get linear response. 

The existence of the adiabatic regime is obviously a 
quantum mechanical effect. If we take the formal limit 
h the adiabaticity condition F ^ A breaks down. 
In fact the proper classical limit is non-perturbative, be- 
cause also the weaker condition F ^ Af, does not sur- 
vive the h ^ limit. For further details see [16-19]. 
In the non-perturbative regime the quantum mechani- 
cal derivation of Kubo formula is not valid. Indeed we 
have demonstrated [18] the failure of Kubo formula in 
case of random-matrix models. But if the system has a 
classical limit, then Kubo formula still holds in the non- 
perturbative regime due to semiclassical (rather than 
quantum- mechanical) reasons. 

The discussion of dissipation assumes a generic situa- 
tion such that the Schrodinger equation does not have 
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a stationary solution. This means that driven one- 
dimensional systems arc automatically excluded. An- 
other non-generic possibility is to consider a special driv- 
ing scheme, such as translation, rotation or dilation [23] . 
In such case the time dependent Hamiltonian TC{x{t)) 
possesses a stationary solution (provided the "velocity" 
X is kept constant). Consequently we do not have a dis- 
sipation effect. In Section 11 we discuss the simplest 
example of pumping by translation, where the station- 
ary adiabatic solution of Section 5 is in fact exact, and 
no dissipation arises. 



IX. APPLICATION TO PUMPING 

So far we have discussed the response for driving in a 
very general way. From now on we focus on a system with 
a ring geometry as described in the Introduction, and 
illustrated in Fig. 1. The shape of the ring is controlled 
by some parameters xi and 2:2, and Xs is the magnetic 
flux. The generalized force which is conjugate to the 
flux is the current. The time integral over the current is 
the transported charge: 



Q 



j{F^)dt 



(38) 



In fact a less misleading terminology is to talk about 
"probability current" and "integrated probability cur- 
rent". From a purely mathematical point of view it is 
not important whether the transported particle has an 
electrical charge. 

Disregarding a possible persistent current contribu- 
tion, the expression for the pumped charge is: 



Q = 



G ■ dx 



BAdx 



(39) 



If we neglect the flrst term, which is associated with the 
dissipation effect, and average the second ("adiabatic") 
term over the flux, then we get 



1 

'2^ 



B • dec A da; 



(40) 



The integration should be taken over a cylinder of ver- 
tical height 2Trh, and whose basis is determined by the 
projection of the pumping cycle onto the {x-i,X2) plane. 

We already pointed out that the Berry phase 
{l/h) § An ■ dx is gauge invariant. Therefore from Stokes 
law it follows that (l/7i) J B ■ dx A dx is independent of 
the surface, and therefore {1/K)§ B ■ dx Adx with closed 
surface should be 27rx integer. Integrating over a cylin- 
der, as in Eq.(40), is effectively like integrating over a 
closed surface (because of the 27r periodicity in the ver- 
tical direction). This means that the flux averaged Q of 
Eq.(40) has to be an integer. 

The common interest is in pumping cycles in the <E> = 
plane. This means that the zero order conservative con- 



tribution to Q, due to a persistent current, does not ex- 
ist. Furthermore, from the reciprocity relations (see Ap- 
pendix B) it follows that G-^i = -G^^ and = -G^^, 
which should be contrasted with G^^ = G^^ . This means 
that a pumping cycle in the $ = plane is purely adia- 
batic: there; is no dissipative contribution to Q. Only the 
B field (second term in Eq.(39)) is relevant to the calcu- 
lation of the pumped charge, and its vertical component 
B^^ vanishes due to the time reversal symmetry. 

The absence of dissipative contribution for a cycle in 
the $ = plane, does not imply that dissipation is not 
an issue. The symmetric part of the conductance ma- 
trix rf^ is in general non-zero, leading to an energy ab- 
sorption rate which is proportional to x^ . This implies 
that the energy absorption per cycle is proportional to 
Therefore we are able to minimize the dissipation 
effect by making the pumping cycle very slow. Further- 
more, if we get into the quantum- mechanical adiabatic 
regime, then rf^ becomes extremely small, and then we 
can neglect the dissipation effect as long as quantum- 
mechanical adiabaticity can be trusted. 

Whenever the dissipation effect cannot be neglected, 
one should specify whether or how a stationary operation 
is achieved. In case of pumping in open system the sta- 
tionary operation is implicitly guaranteed by having equi- 
librated reservoirs, where the extra energy is dissipated 
to infinity. In case of pumping in closed system the issue 
of stationary operation is more subtle: In the adiabatic 
regime, to the extend that adiabaticity can be trusted, 
we have a stationary solution to the transport problem, 
as defined in Section 5. But outside of the adiabatic 
regime wc have diffusion in energy space (Appendix E) 
leading to a slow energy absorption (dissipation). Thus 
a driven system is heated up gradually (though possibly 
very slowly). Strictly speaking a stationary operation 
is not achieved, unless the system is in (weak) thermal 
contact with some large bath. Another way to reach a 
stationary operation, that does not involve an external 
bath, is by having an effectively bounded phase space. 
This is the case with the mixed phase space example 
which is discussed in Ref.[12]. There the stochastic-like 
motion takes place in a bounded chaotic region in phase 
space. 



X. CLASSICAL DISSIPATIVE PUMPING 

Before we discuss the quantum mechanical pumping, 
it is instructive to bring simple examples for classical 
pumping. In the following we consider one particle (r) 
in a two dimensional ring as in Fig. la. 

The first example is for classical dissipative pumping. 
The conductance G = G^^ can be calculated for this 
system [24] leading to a mesoscopic variation of the Drudc 
formula. The current is given by Ohm law / = — G x $, 
where — $ is the electro- motive- force. 

Consider now the following pumping cycle: Change 
the fiux from $1 to $2) hence pumping charge 
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Q = —G{1) X (<I>2 — $i). Change the conductance from 
G{1) to G(2) by modifying the shape of the ring. 
Change the flux from <I>2 back to $i, hence piuiip- 
ing charge (5(2) = -G(2) x ($i - $2)- Consequently the 
net pumping is 



XII. QUANTUM PUMPING 

We turn now to the quantum mechanical case. Con- 
sider an adiabatic cycle that involves a particular energy 
level n. This level is assumed to have a degeneracy point 



Q = (G(2) - G(l)) X ($2 - *i) 



(41) 



at {x\ 



(0) ^(0) 



, $(°)). It follows that in fact there is a ver- 



tical chain of degeneracy points: 



Thus we have used the dissipativc part of the conduc- 
tance matrix (first term in Eq.(39)) in order to pump 
charge. In the quantum mechanical version of this exam- 
ple extra care should be taken with respect to the zero 
order contribution of the persistent current. 



XI. CLASSICAL ADIABATIC PUMPING 

The second example is for classical adiabatic pumping. 
The idea is to trap the particle inside the ring by a poten- 
tial well, and then to make a translation of the trap along 
a circle. The result of such a cycle is evidently Q = 1. 
We would like to see how this trivial result emerges form 
the Kubo formula. 

Let (r,p) be the canonical coordinate of the particle 
in the ring, while {xi,X2) are the center coordinate of a 
trapping potential. The Hamiltonian is: 



chain 



(0) „(0) 

tJbO 



) ^ $(0) + 27r?i X integer) (46) 



n{r,p;x{t)) = 



1 



pI + 



P\\ 



2n^x'f + x\ 



2m 

+ U,r^p{ri-xi{t),r2-X2{t)) 



-Mt) 



(42) 



where p\\ and p± are the components of the momentum 

along the ring and in the perpendicular (transverse) di- 
rections. The pumping is done simply by cycling the po- 
sition of the trap. The translation of the trap is assumed 
to be along an inside circle of radius R, 



These degeneracy points are important for the geo- 
metrical understanding of the B field, as implied by 
Eq.(18). Every degeneracy point is like a monopole 
charge. The total flux that emerges from each monopole 
must be 27r?ix integer for a reason that was explained af- 
ter Eq.(40). Thus the monopoles are quantized in units 
of Ji/2. The B field which is created (so to say) by a ver- 
tical chain of monopoles may have a different near field 
and far field behavior, which we discuss below. 

The far field region exists if the chains are well isolated. 
Later we explain that "far" means gx ^ 1, where gx is 
the Thouless conductance. The far field is obtained by 
regarding the chain as a smooth line. This leads qualita- 
tively to the same field as in Eq.(45). Consequently, for 
a "large radius" pumping cycle in the $ = plane, we 
get IQI « 1. In the following we are interested in the de- 
viation from exact quantization: If 0^"' = we expect to 
have IQI > 1, while if 0^°' = tt we expect \Q\ < 1. Only 
for the (j) averaged Q of Eq. (40) we get exact quantization. 

The deviation from \Q\ « 1 is extremely large if we 
consider a tight pumping cycle around a 0^°^ = degen- 
eracy. After linear transformation of the shape param- 
eters, the energy splitting A = — E^n of the energy 
level n from its neighboring (nearly degenerated) level m 
can be written as 



A = {{x^-xf^f + (X2-4°T + c'{^-cf>^'^r)y' (47) 
x{t) = {Rcos{flt), Rsm{flt),^=const) (43) where c is a constant. The monopole field is accordingly 



In this problem the stationary solution of Section 5 is 
an exact solution. Namely 



B = ±- 



(xi 



.(0) 



(0) 

,X2-X), ',X3-Xl ') 



(0)^ 



\n{xm 



(44) 



where \n{x)) 1-^ tp^^'>{r — x) are the eigenfunctions of a 
particle in the trap. Eq.(44) is just Galilei transformation 
from the moving (trap) frame to the Laboratory frame. 

It is a-priori clear that in this problem the pumped 
charge per cycle is Q = 1) irrespective of Therefore 
the B field must be 



2 {{x^-xf^r + {x2-xrr + (f )'(a;3-4 T)'/ 



where the prefactor is determined by the requirement 
of having a single {fi/2) monopole charge. Assuming a 
pumping cycle of radius R in the $ = plane we get 
from the second term of Eq.(39) 



#8) 



Q 



B 



(xi, 2:2,0) 



where 



2tt{xI 



4) 



(45) 



B Arfa; 



1 d^A 



9t 



I ■ : I - can be verified by calculation via Eq.(18). The sin- 
gularity along the 0:3 axis is not of quantum mechanical 
origin: It is not due to degeneracies, but rather due to 
the diverging current operator {dH/dxs ex 1/ ^/xl+^). 



R^ 



(49) 



(50) 



is a practical definition for the Thouless conductance in 

this context. It is used here simply as a measure for the 
sensitivity of an energy level to the magnetic fiux 
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What we want to do in the next sections is to interpo- 
late between the near field result, which is Q = 0{yfgT), 
and the far field result, which is Q = 0(1). For this pur- 
pose it is convenient to consider a particular model that 
can be solved exactly. 

XIII. THE DOUBLE BARRIER MODEL 

A simple example for quantum pumping is the double 
barrier model. An open geometry version of this model 
has been analyzed in [8] using the S matrix approach. 
The analogous closed geometry version is obtained by 
considering a one-dimensional ring with two delta bar- 
riers. As we arc going to explain below, the pumping 
process in this model can be regarded as a particular ex- 
ample of an adiabatic transfer scheme: The electrons are 
adiabatically transfcred from state to state one by one as 
in "musical chair game" . 

The two delta barriers version of the double barrier 
model is illustrated in Fig. 2. The length of the ring is 
L, with periodic boundary conditions on — (L/2) < r < 
{L/2). A dot region \Q\ < a/2 is defined by the potential 

C/(.;c„c.) = l^(r + |)+l^(.-^) (51) 

It is assumed that ci and C2 are small enough so one 
can classify the ring eigenstates into two categories: wire 
states, and dot states. The latter arc those states that 
are localized in the dot region \Q\ < a/2 in the limit of 
infinitely high barriers. We define the Fermi energy as 
the energy of the last occupied wire level in the limit of 
infinitely high barriers. 

The three parameters that we can control are the flux 
2:3 = $ = h(j), the bias xi = ci — C2, and the dot potential 
X2 = i?dot which is related to ci -f C2. The energy E^^-i 
correspond to the dot state which is closest to the Fermi 
energy Ep from above. We assume that the other dot 
levels are much further away from the Fermi energy, and 
can be ignored. Note that another possible way to control 
the dot potential, is simply by changing a gate voltage: 
That means to assume that there is a control over the 
potential floor in the region \Q\ < a/2. 

The pumping cycle is assumed to be in the $ = plane, 
so there is no issue of conservative persistent current con- 
tribution. The pumping cycle is defined as follows: We 
start with a positive bias {xi > 0) and lower the dot po- 
tential from a large X2 > Ep value to a small X2 < Ep 
value. As a result, one electron is transfered via the 
left barrier into the dot region. Then we invert the bias 
(a;i < 0) and raise back X2- As a result the electron is 
transfered back into the wire via the right barrier. 

A closer look at the above scenario (Fig. 2b) reveals the 
following: As we lower the dot potential across a wire 
level, an electron is adiabatically transfered once from 
left to right and then from right to left. As long as the 
bias is positive (a;i > 0) the net charge being pumped 
is very small {\Q\ <C 1). Only the lowest wire level that 



participate in the pumping cycle carries Q = 0(1) net 
charge: It takes an electron from the left side, and after 
the bias reversal it emits it into the right side. Thus 
the pumping process in this model can be regarded as a 
particular example [14] of an adiabatic transfer scheme: 
The electrons are adiabatically transfered from state to 
state, one by one, as in "musical chair game" . 

For a single occupied level the net Q is the sum of 
charge transfer events that take place in four avoided 
crossings (two avoided crossings in case of the lowest 
level). For many particle occupation the total Q is the 
sum over the net Qs which arc carried by individual lev- 
els. For a dense zero temperature Fermi occupation the 
summation over all the net Qs is a telescopic sum, leav- 
ing non-canceling contributions only from the first and 
the last adiabatic crossings. The latter involve the last 
occupied level at the Fermi energy. 



XIV. THE THREE SITE LATTICE 
HAMILTONIAN 

Rather than analyzing the two-delta-barriers version 
of the double barrier model, we consider below a simpli- 
fied version that still contains the same essential ingredi- 
ents. This is obtained by considering a three site lattice 
Hamiltonian. The advantage is obviously the possibil- 
ity to make an exact analytical treatment that does not 
involve approximations. 

The middle site in the three site lattice Hamiltonian 
supports a single dot state, while the two other sites sup- 
port two wire states. The Hamiltonian is 

/ Cl e'-^X 
H Cl u C2 (52) 

Ve-'"^ C2 / 

The three parameters are the bias Xi = Ci — C2, the dot 
energy X2 = u, and the flux X3 = $ = hcj). For presenta- 
tion purpose we assume that < Ci, C2 -C 1, and charac- 
terize the wire-dot coupling by the parameter c = ,Jc\C2. 

The eigenstates are En- Disregarding the interaction 
with the dot (c = 0) we have two wire states with 
E = ±1. This implies degeneracies for 0:2 = u = =f1- 
Once we switch on the coupling (c > 0), the only possible 
degeneracies are between the even dot state and the odd 
wire state of the mirror symmetric Hamiltonian {x\ = 0). 
The flux should be either integer (for degeneracy of the 
dot level with the lower wire level), or half integer (for 
degeneracy of the dot level with the upper wire level). 
Thus we have two vertical chains of degeneracies: 

The negative chain = (0, — l-|-c^, 27r7i x integer) 
The positive chain = (0, -|-1— c^, tt -|- 27r?i x integer) 

In order to calculate the B field and pumped charge 

Q, we have to find the eigenvalues and the eigenvectors 
of the Hamiltonian matrix. The secular equation for the 
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eigenvalues is 

- uE'^ - (1 + cf + (^)E + u- 2ciC2 cos(^) = 
Using the notations 



27 



cos(6') = 



i (1 + Ci + c^u - + C1C2 cos{(j)) 

n 



the roots of the above cubic equation are: 

1 



E„ 



2^5 cos ( ^^' + "Y 



(53) 



where n = 0, ±1. The corresponding eigenstates are: 



1 



v-^ Vcie-*'^ + C2i;„ 



where 5 is the normalization, namely 

S = {l-Elf + {c^+C2E^f + {c2+c^E^f 



(54) 



(55) 



For the calculation of the pumped charge in the next 
paragraph it is useful to notice that for i? = ±1 the 
normalization is 5 = 2(ci ± 02)^, while for = the 
normalization is 5 ~ 1. 

After some algebra wc find that the first component of 
the B field in the $ = plane is 



\ ou 



—n[x) 



-{cl-cD 



(56) 
(57) 



Which is illustrated in Fig. 3. From here it follow that if 
we keep constant bias, and change only X2 = u, then the 
pumped charge is: 



Q 



= - j B'dx2=-{cl-cl)^ 



(58) 



For a planar ($ = 0) pumping cycle around the negative 
vertical chain the main contribution to Q comes from the 
two crossings of the a;2 ~ — 1 line. Hence we get 



Q 



Cl + C2 
Cl - C2 



(59) 



where the Thoulcss conductance in this context refers to 
the avoided crossing, and is defined as 



9t 



1 a^A 



4>=o 



2ciC2 
(Cl - C2f 



(60) 



A similar calculation of the pumped charge for a planar 
cycle around the positive chain leads to 



Q 



Cl - C2 
Cl + C2 



VI -257 



(61) 



with qt = 2ciC2/(ci + 02)^. In both cases we have ap- 
proximate quantization Q = ±1 + 0{gT) for ^ 1, 
while for a tight cycle cither Q ^ 00 or Q ^ de- 
pending on which line of degeneracies is being encircled. 
If the pumping cycle encircles both chains then we get 
Q = 4ciC2/(cf — c|). In the latter case Q = 0{gT) for 
5t "C 1, with no indication for quantization. 



XV. SUMMARY AND DISCUSSION 

Wc have shown how the Kubo formalism can be used 
in order to derive both classical and quantum mechanical 
results for the pumped charge Q in a closed system. In 
this formulation the distinction between dissipative and 
non-dissipative contributions is manifest. 

Within the framework of the Kubo formalism (disre- 
garding non-linear corrections) we have made a distinc- 
tion between the following levels of treatment: 

• Strict adiabaticity 

(outcome of zero order treatment) 

• Adiabatic transport 

(outcome of stationary first order treatment) 

• Dissipation 

(the result of first order transitions) 

In the adiabatic regime one can assume a stationary so- 
lution to the adiabatic equation, which implies no dis- 
sipation effect. This leads to the picture of adiabatic 
transport, where the Berry phase is the outcome of a 
zero order treatment, while the "geometric magnetism" 
of Eq.(26) is the outcome of a first order treatment of the 
inter-level couplings. 

In some very special cases (translations, rotations and 
dilations) this assumption (of having a stationary solu- 
tion) is in fact exact, but in generic circumstance this 
assumption is an approximation. Outside of the adia- 
batic regime the stationary solution cannot be trusted. 

Assuming quantized chaotic dynamics one argues 
that Fermi-golden-rule transitions between levels lead to 
(slow) difi^usion in energy (Eq.(El)). This leads to the 
emergence of the dissipative part in the Kubo formula. 
We have obtained an expression (Eq.(35)) for the en- 
ergy scale F oc l.-Ep/'"' that controls the dissipative ef- 
fect. We have explained that the dissipative contribu- 
tion to the Kubo formula is valid only in the regime 
A < r <C Afc. Otherwise the dynamics is either of adia- 
batic nature (F <C A) or non-perturbative (F > A). 

In order to calculate the pumped charge Q we have 
to perform a closed line integral over the conductance 
(Eq.(39)). This may have in general both adiabatic and 
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dissipative contributions. For the common pumping cy- 
cle in the $ = plane, only the adiabatic contribution 
exists. This follows from the reciprocity relations (Sec- 
tion 9). Still we have emphasized (without any contradic- 
tion) that in the same circumstances a dissipation effect 
typically accompanies the pumping process. 

The quantum adiabatic contribution to the pumping 
is determined by a line integral over a B field which is 
created by m,onopoles. The monopolcs, which are related 
to the degeneracies of the Hamiltonian, are located along 
vertical chains in x space (Eq.(46)). The 3 site model 
provides the simplest example for such vertical chains: 
By calculating the B field which is created (so to say) by 
these chains, we were able to determine the charge which 
is pumped during a cycle (e.g. Eq.(59)). 

The (monopoles of the) vertical chains have near field 
regions (Eq.(48)). If the chains are well isolated in x 
space, then there are also far field regions. The far field 
regions are defined as those where the Thouless conduc- 
tance is very small {gr ^ !)• Pumping cycles that are 
contained in the far field region of a given chain lead to an 
approximately quantized pumping Q = integer -I- 0(gT)- 
It is important to realize that the existence of far field 
regions in x space is associated with having a low di- 
mensional system far away from the classical limit. In a 
quantized chaotic system it is unlikely to have <C 1 
along a pumping cycle. As we take the ?l ^ limit the 
vertical chains become very dense, and the far field re- 
gions disappear. 

In the subtle limiting case of open geometry we expect 
to get agreement with the S'-matrix formula of Biittiker 
Pretre and Thomas (BPT) [6]. Using the notations of 
the present Paper the BPT formula for the current that 
comes out of (say) the right lead can be written as: 

= ^trace ( P^S'^\ (62) 
2m \ dxj J ^ ' 

where P is the projector on the right lead channels. For 
G^^ the above reduces to the Landauer formula. The 
details regarding the relation between the Kubo formula 
and the BPT formula will be published in a separate 
paper [10]. Here wc just note that the derivation is based 
on a generalization of the Fisher-Lee approach [3, 25, 26]. 

Finally it is important to remember that the theory 
of driven systems is the corner stone for the analysis of 
interaction between "slow" and "fast" degrees of freedom. 
Assume that that the x^ are in fact dynamical variable, 
and that the conjugate momenta are pj. The standard 
textbook example is the study of diatomic molecules. In 
such case x^ are the locations of the nuclei. The total 
Hamiltonian is assumed to be of the general form 

W,„,., = ^^p2 ^ H(a;) (63) 

3 

where W is the Hamiltonian of the "fast" degrees of free- 
dom (in the context of molecular physics these are the 
electrons). Rather than using the standard basis, one can 



use the Born-Oppenheimer basis \x,n(x)) = \x)(^\n{x)). 
Then the Hamiltonian can be written as 

3 

where the interaction term is consistent with Eq.(19). 
Thus it is evident that the theory of driven systems is 
a special limit of this problem, which is obtained if we 
treat the classical variables. 
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APPENDIX A: THE KUBO FORMULA: 
STANDARD DERIVATION 



APPENDIX B: REMARKS REGARDING THE 
GENERALIZED SUSCEPTIBILITY 



In this Appendix wc present an elementary textbook- 
style derivation of the Kubo formula. For notational sim- 
plicity we write the Hamiltonian as 7i = Hq — /(t)V. It 
is assumed that the system, in the absence of driving, 
is prepared in a stationary state po- In the presence of 
driving wc look for a first order solution p{t) = po + p{t)- 
The equation for p{t) is: 



dm 

dt 



-z[Ho,p(t)] + i/(t)[V,po] 



(Al) 



This equation can be re-written as 

-(Uo(i)-V(i)Uo(t)) « i/(t)[Uo(t)-^VUoW,po] 

where Uo{t) is the evolution operator which is generated 
by Hq. The solution of the latter equation is 



m 



j\[M{-{t-t')),p^]f{t')dt' (A2) 



where we use the usual definition of the "interaction pic- 
ture" operator V(t) = Uo(t)-i VUo(t). 

Consider now the time dependence of the expectation 
value (F)t = trace(Fp(t)) of an observable. Disregarding 
the zero order contribution, the first order expression is 

{?)t « y*ztrace(F[V(-(t-t')),/Oo]) f{t')dt' 
a{t - 1') f{t')dt' 



-I 



where the response kernel a(r) is defined for r > as 

a(r) = z trace (F[\/(—r), po]) 

= i trace ([F,V(-t)]/9o) 

= i([F,V(-r)]) 

= i([F(r),V]) (A3) 

We have used above the cyclic property of the trace oper- 
ation; the stationarity UoPoUq ^ = Po of the unperturbed 
state; and the definition F(t) = Uo(t)^^FUo(t). 



In this appendix we would like to further illuminate 
the relation between the generalized susceptibility and 
the conductance matrix. The generalized susceptibility 
x'^-^(a;) is the Fourier transform of the causal response 
kernel a''^ {t). Therefore it is an analytic function in the 
upper half of the complex oj plan, whoso real and imag- 
inary parts are related by Hilbert transforms (Kramers- 
Kronig relations): 



-< 



Im[x'=^(w')] du' 



(Bl) 



The imaginary part of x'^-' (w) is the sine transforms of 
Q!*^-' (t), and therefore it is proportional to u for small 
frequencies. Consequently it is convenient to write the 
Fourier transformed version of Eq.(3) as 

[{F')U = Y^Xo^i^n^jh-f^^'i^nijU (B2) 



where the dissipation coefficient is defined as 

/.(.;) = = r a^^irf^^dr (B3) 

U Jo U) 

In this paper we ignore the first term in Eq.(B2) which 
signify the non-dissipative in-phase response. Rather 
we put the emphasis on the "DC limit" {lo 0) 
of the second term. Thus the conductance matrix 
Qkj _ ^kj^^^ _^ jg j^g^ synonym for the term "dis- 
sipation coefficient". However, "conductance" is a bet- 
ter (less misleading) terminology: it does not have the 
(wrong) connotation of being specifically associated with 
dissipation, and consequently it is less confusing to say 
that it contains a (non-dissipative) adiabatic component. 

For systems where time reversal symmetry is broken 
due to the presence of a magnetic field B, the response 
kernel, and consequently the generalized susceptibility 
and the conductance matrix satisfies the Onsager reci- 
procity relations 



a'^(r,-S) 



[±] G^'iB) 



(B4) 
(B5) 
(B6) 



where the plus (minus) applies if the signs of and 
transform (not) in the same way under time reversal. 
These reciprocity relations follow from the Kubo formula 
(Eq.(28), using K'i{-T,-B) = -[±]K'^ (t, B), together 
with the trivial identity K^^{-t,B) = -K^^{t,B). In 
Section. 9 we discuss the implications of the reciprocity 
relations in the context of pumping. 
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APPENDIX C: EXPRESSIONS FOR B AND G 

The functions C"-' (r) and K'^^ (r) are the expectation values of hermitian operators. Therefore they are real functions. 
It follows that the real part of their Fourier transform is a symmetric function with respect to cj, while the imaginary 
part of their Fourier transform is anti symmetric with respect to w. By definition they satisfy C*-' (r) = C^^{~t) and 
K^^{t) = —K^^{—t). It is convenient to regard them as the real and imaginary parts of one complex function $'-' (t). 
Namely, 



$'J'(t) = {F\t)F^{0)) 
C'\t) = i($'^(T) + $^^(-r)) 
K'\t) = |($'^(r)-$^'(-T)) 



C'^{T)-i-^K'^{T) 



(CI) 
(C2) 
(C3) 



It is possible to express the decomposition G*-' = r;'^ + in terms of K'^^{uj). Using the definition Eq.(29) we get: 



Jo 



K'^{T)TdT 



--f 



llm[if*^(w)] 
2 w 



(C4) 



01=0 



The first term is antisymmetric with respect to its indexes, and is identified as i?*-' . The second term is symmetric 
with respect to its indexes, and is identified as if^ . The last step in the above derivation involves the following identity 
that hold for any real function /(r) 



2-K 



1 Ini[/(a>)] 

2 UJ 



(C5) 



w=0 



Note that Im[/(a;)] is the sine transform of /(r), and therefore it is proportional to w is the limit of small frequencies. 

It is of practical value to re-derive Eq.(C4) by writing $'-' (r) using the energies En and the matrix elements F^^. 
Then we can get from it straightforwardly (using the definitions) all the other expressions. Namely: 



$y(r) 
l>'^(a;) 



n m ^ 

J2fiEn)J2F^^Fin 2nd (cj 

n 



h 



n 



nm mn 



+ 



P3 pi 
nm mn 



^%uj-{E^-En)+iQ nLJ+{Em-En)+iO 

-2TThJ2fiEn) J2 ^[FnmFlnl S' {Em ' E^) 



2nY^f{E^) E ^^^^^ 

i-L \ \^m ^nj 
n m(^n) 



(C6) 

(C7) 

(C8) 
(C9) 

(CIO) 



One observes that the expression for i?*-' coincides with the adiabatic transport result Eq.(26). Alternatively this 
identification can be obtained by expressing the sum in Eq.(ClO) as an integral, getting form it the first term in 
Eq.(C4): 



2 r 



Im[^'3{oj)] du> 



2 

n. 



lm[C'3{iv)]-^Re[K'i{u)]duj 



27r 



2^ 



(Cll) 
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APPENDIX D: EXPRESSING K{uj) USING C{u!) 
We can use the following manipulation in order to relate ^'■' (w) to C"^ (a;), 

K'\u) = ^/(i?„)i^;^H (Di) 

n 

nm 



nm 



nm 



nni 



n 

where we use the notation Wnm = {En — Em)/h. The third line differs from the second line by permutation of the 
dummy summation indexes, while the fourth line is the sum of the second and the third lines divided by 2. In the last 
equality we assume small w. If the levels are very dense, then we can replace the summation by integration, leading 
to the relation: 

J g{E)dE f{E) k%{uj) = -iu J g{E)dE f{E) C%{u) (D2) 

where K^^{ui) and C'^(w) are microcanonically smoothed functions. Since this equality hold for any smoothed f{E), 
it follows that the following relation holds (in the limit w — > 0): 

K'ii.:) = [giEKiico)] (D3) 

If we do not assume small w, but instead assume canonical state, then a variation on the last steps in Eq.(Dl), using 
the fact that {f{En)—f{Em,))/{f{En)+f{Em)) = tanh((i?„— £'m)/(2T)) is an odd function, leads to the relation 

k'^{u) = iwxr^tanhf^) C'^iu) (D4) 



Jiu \2T J 

Upon substitution of the above expressions in the Kubo formula for t?'-' , one obtains the Fluctuation-Dissipation 
relation. 



APPENDIX E: THE KUBO FORMULA AND THE 
DIFFUSION IN ENERGY SPACE 

The illuminating derivation of Eq.(33) is based on the 
observation that energy absorption is related to having 
diffusion in energy space. Let us assume that the proba- 
bility distribution p{E) = g{E)f{E) of the energy satis- 
fies the following diffusion equation: 



The energy of the system is {H) = J Ep{E)dE. It follows 
that the rate of energy absorption is 

For a microcanonical preparation we get 

This diffusion-dissipation relation reduces immediately to 
the fluctuation-dissipation relation if we assume that the 
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diffusion in energy space due to the driving is given by 

= iT^^Ei^^^) i^ij (E4) 



Thus it is clear that a theory for Unear response should 
establish that there is a diffusion process in energy space 
due to the driving, and that the diffusion coefficient 
is given by Eq.(E4). More importantly, this approach 
also allows treating cases where the expression for De is 
non-perturbative, while the diffusion-dissipation relation 
Eq.(E3) still holds! 

A full exposition (and further reference) for this route 
of derivation can be found in [16-19]. Here we shall 
give just the classical derivation of Eq.(E4), which is ex- 
tremely simple. We start with the identity 



(E5) 



Assuming (for presentation purpose) that the rates Xk 
are constant numbers, it follows that energy changes are 
related to the fluctuating F'^(t) as follows: 



SE 



{n)t - (H)o 



E 



F\t')dt' (E6) 



Squaring this expression, and performing microcanonical 
averaging over initial conditions we obtain: 



JO 



c'iit" 



t')dt'dt" (E7) 



where C|(i" - t') = {F'{t')F^ {t")) is the correlation 
function. For very short times this equation implies 

"ballistic" spreading (SE'^ cx t^) while on iriterniediate 
time scales it leads to diffusive spreading 5E'^{t) = 2DEt, 
where 



1 

De = ^ l^XiXj 



&i{T)dT 



(E8) 



The latter result assumes a short correlation time. This 
is also the reason that the integration over r can be ex- 
tended form —00 to -|-oo. Hence we get Eq.(E4). We note 
that for long times the systems deviates significantly from 
the initial microcanonical preparation. Hence, for long 
times, one should justify the use of the diffusion equa- 
tion (El). This leads to the classical slowness condition 
which is discussed in Ref.[19]. 
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FIGl. Illustration of a ring system (a). The shape of the ring is controlled by some parameters xi and X2. The flux through the ring 
is = A system with equivalent topology, and abstraction of the model arc presented in (b) and (c). The "dot" can be represented 
by an S matrix that depends on xi and X2- In (d) also the flux xs is regarded as a parameter of the dot. If we cut the wire in (d) we get 
the open two lead geometry of (e). If we put many such units in series we get the period system in (f). 




position tifne 

FIG2. Schematic illustration of quantum pumping in a closed wire-dot system. The net charge via the third level (thick solid line on 
the right) is vanishingly small: As the dot potential is lowered an electron is taken from the left side (first avoided crossing), and then 
emitted back to the left side (second avoided crossing). Assuming that the bias is inverted before the dot potential is raised back, only 
the second level carry a net charge Q = 0(1). 




FIGS. The first component of the B field for a particle in the middle level of the 3 site lattice model. It is plotted as a function of the 

dot potential X2 = u. The other parameters arc = 0, and ci = 0.1, while C2 = 0.04 for the thick line and C2 = 0.02 for the thin line. In 
the limit C2 0, all the charge that is transfered from the left side into the dot during the first avoided crossing, is emitted back into the 
left side during the second avoided crossing. Inset: The eigenenergies E„{x) for the C2 = 0.04 calculation. 



